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Abstract. Previous work by Rubinstein IIRubll and Gao [Gaol computed the n-level den- 
sities for families of quadratic Dirichlet i-functions for test functions 0i , . . . , supported 
^J[) ill X^iLi < 2, and showed agreement with random matrix theory predictions in this 

^ . range for n < 3 but only in a restricted range for larger n. We extend these results and 

show agreement for n < 7, and reduce higher n to a Fourier transform identity. The proof 
involves adopting a new combinatorial perspective to convert all terms to a canonical form, 
which facilitates the comparison of the two sides. 
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1. Introduction 

Assuming the Generalized Riemann Hypothesis (GRH), the non-trivial zeros of L-functions 
lie on the line = 1/2. The distribution of these zeros govern the behavior of a variety of 
problems, ranging from the distribution of primes in arithmetic progressions to the size of 
the class number to the geometric rank of the Mordell-Weil group of elliptic curves, among 
others [ICTl IDai IGoi IGZl IRubSall . In many instances we need to know more than just the 
fact that the zeros lie on the line, but additionally how they are distributed on the line. 

One of the most successful approaches to modeling these zeros is through Random Ma- 
trix Theory. Originally arising in statistical investigations HWisH . the subject flourished in 



the 1950s and 1960s with the work of Wigner [Wig 1] [Wipl |Wig3[ [Wig4l [Wig5 ] , Dyson 



[ |Dyl[ |Dy2[ and others, who applied it to describe the energy levels of heavy nuclei. In 
the 1970s, Montgomery and Dyson HMonl noticed that the 2-level correlation of zeros of 
the Riemann zeta function matched those of the Gaussian Unitary Ensemble (GUE); see 
HHai IFirM i I for more on the history. Since then Random Matrix Theory has made precise 
statements about the main term in the behavior of numerous statistics involving zeros of 
L-functions llConl IKeSnli IKeSn2[ IKeSn3ll . 

While the limiting behavior of n-level correlations of a single L-function have been 
shown to agree (for suitable test functions) with the scaling limit of the GUE [ jHqj |Mon[ 
IRSL the behavior near the central point is different for different L-functions, and depends 
on the arithmetic of the form (for example, the order of vanishing of L-functions attached 
to elliptic curves is conjecturally equal to the rank of the Mordell-Weil group). To study 
these low-lying zeros, Katz and Sarnak UKaSall IKaSa21l introduced the n-level density. 
Assuming GRH, the non-trivial zeros of an L-function L(s, /) are 1/2 + i'y^^^^ with 7^'' 

real, where ■ ■ ■ < 7j < 7/"^'' < 7/^'' < if^ < ■ ■ ■ if the sign of the functional equation 

is even (if it is odd, there is an extra zero: 7^-°^ = 0). The n-level density for a finite family 
of L-functions is 

^'"i^: ^) - 4 E E O^^r) • • ■ C^r) . (>■» 

where the (pi are even Schwartz functions whose Fourier transforms have compact sup- 
port and logi? is a normalization parameter (essentially the average of the logarithms of 
the analytic conductors) so that the scaled zeros near the central point have mean spacing 
1. The Katz-Samak Density Conjecture states that as the conductors tend to infinity the 
distribution of the scaled zeros near the central point converges to the same limiting distri- 
bution as the normalized eigenvalues near 1 of a subgroup of the unitary group U (N) as 
— 7- 00. The corresponding group is typically unitary, symplectic, or orthogonal matrices 
(or a trivial modification to take into account forced zeros at the central point). There is 
strong evidence for this conjecture. First, in the function field case the correspondence is 
clear as the subgroup is the monodromy group. Second, there are now many families of 
L-functions where we can prove agreement for suitably restricted test functions, includ- 
ing Dirichlet L-functions, elliptic curves, cuspidal newforms, Maass forms, number field 
L-functions, symmetric powers of GL2 automorphic representations and Rankin-Selberg 
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convolutions of the above, to name a few IIAILMZi iDMTl lDM2l IFioMii ED iGaol iGUl iHMl 
IHmiimiKaSnilKaSa2imiMllPel[OSTllOS2lim 

In this paper we study the low-lying zeros (i.e., those near the central point) of quadratic 
Dirichlet L-functions via the n-level density. In his thesis Rubinstein [Rub J showed these 
agree with the scaling limit of symplectic matrices whenever 0i, . . . , are supported in 
Yll=i < 1- Gao UGaol extended this result in his thesis. It is important to have as large 
support as possible, as frequently extending the support is related to finer questions about 
the arithmetic of the family. Interestingly, while Gao was able to compute the number the- 
ory side for test functions supported in Yll=i Wi\ < 2, he was only able to show agreement 
with the Katz-Samak determinantal expansion for the symplectic ensemble forn < 3. 

This created an annoying situation in the literature, where both number theory and ran- 
dom matrix theory had been computed in the regime Yll=i l^il < 2, but could only be 
shown to agree in this full range for n < 3. Gao's proof involved using ad hoc Fourier 
transform identities to match the manageable number of terms present for such small n. 
Unfortunately, the number of summands grows very rapidly with n, and this approach be- 
comes impractical for higher n. 

In this paper, we further extend the agreement between number theory and random ma- 
trix theory. Our proof is in two steps. First, we resolve a combinatorial obstruction by 
rewriting both densities using the same combinatorial perspective: we express the terms of 
the densities in terms of certain pairs of set partitions. This allows us to show agreement 
between most of the terms arising in the densities, for any n. Second, we reduce the Den- 
sity Conjecture (in the range Yll.=i Wi\ < 2) to showing that a term arising in the random 
matrix theory is the Mobius transform (over the lattice of set partitions) of a corresponding 
term from number theory. We cannot prove this identity for all n, but we use Mobius in- 
version and properties of Fourier transforms to give it a canonical form that is possible to 
check with a computation. As an application, we verify it for n < 7: 

Theorem 1.1. Let (pi, . . . , (pn be even Schwartz functions with the (pi supported in Yl^=i l^il 
< 2. For n < 7, the Katz-Sarnak Density Conjecture holds for the low-lying zeros of 
quadratic Dirichlet L-functions {L{s, X8d)}> where d is odd and square-free. 

In the above theorem, following Gao UGaol we restricted the family of quadratic char- 
acters. This simplifies the analysis by excluding X2, and facilitates applications of Poisson 
summation in Gao's thesis UGaoL Note that xsd is a real primitive character with even sign 
(i.e., X8d(-1) = !)• 

We briefly sketch the proof. Both sides are known for Yll=i 1^*1 < 2 by URubi IGaol ; the 
difficulty is showing that the two expressions are equal. We proceed as follows. 

(1) We regroup the terms in the random matrix theory in terms of pairs of set partitions 
F, G, such that F_ refines G_ and each block of G is a union of at most two blocks 
of F. 

(2) We do the same to the number theory; this step is more involved because the count- 
ing is naturally 'backwards' there, so the main step is to switch the order in which 
the pairs of partitions are counted. 
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(3) We separate the remaining non-matching terms from the rest of the sum, and show 
that they are all instances of a single Mobius inversion identity. 

(4) We (Mobius-)invert the identity and use properties of Fourier transforms to convert 
all the terms to integrals over ]R>q. We reduce to showing that the integrands are 
identically equal in the region Ui > 0, J2^=i < 2. 

(5) We reduce to a formal polynomial over the subsets of {1, . . . , n}, modulo two re- 
lations that encode the support restriction: this gives an algorithm for showing the 
Fourier identity, which we use to verify up to n = 7. 

Remark 1.2. This work is an extension of the first-named author's 2011 senior thesis at 
Williams College ULevll . There agreement was shown for n < 6 through a more compu- 
tational approach. In the course of extending these results and preparing this manuscript, 
we learned of the work of Entin, Roddity-Gershon and Rudnick HER-GRI . who are able 
to show agreement for all n. Instead of taking a combinatorial approach, they proceed by 
going through a function field analogue and using the limit of large finite fields where the 
hyperelliptic ensemble is shown to have USp statistics. In particular, their results imply 
that our identity holds for all n; it would be interesting to complete the ideas of this paper 
and derive a purely combinatorial proof of this fact. 

The paper is organized as follows. We assume the reader is familiar with URubllGaoll . and 
we will just quote the number theory and random matrix expansions from these works. In 
Owe review some notation and derive some combinatorial results which allow us to recast 
our problem as a related Fourier transform identity. We briefly discuss the obstruction 
which restricts our theorem ton < 7, and see why the two sides at first look so different. We 
continue in the next section by recasting the random matrix and number theory expansions 
to a more amenable form, reducing the problem to the aforementioned Fourier transform 
identity, which we analyze in HI There we rewrite everything in a more tractable canonical 
form, and discuss the verification for n < 7, which completes the proof of Theorem 1 1.1[ 

2. Combinatorial preliminaries 

The purpose of this section is to set the notation for the subsequent combinatorial analy- 
sis, and highlight the technical issues. 

2.1. Set partitions. We recall some basic properties of set partitions. A partition F_ of 
a finite set S* is a collection of subsets F_ = {Fi, . . . , F^} C V{S), such that the Fj are 
nonempty and pairwise disjoint, and S = IJi=i ^i- The Fj are called the blocks of F_ and 
the number k = ^{F) of blocks is the length of F. The set of all partitions of a set S is 
denoted n(S'); when n G N, by abuse of notation we write Y[{n) for n({l, . . . , n}). 

We partially-order 11 (S") by partition refinement: F ^ G if each block of F is contained 
in some block of G_ (equivalently, each block of G is a union of blocks of F). We write 
0_ = {{1}, . . . , {n}} and iV = {{1, . . . , n}} for the minimal and maximal partitions. 

We associate to any partially ordered set P the incidence algebra 



A = {f : P X P ^ C \ f{x, y) = unless x ^ y}, 



(2.1) 
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with pointwise addition and multiplication defined by the convolution *: 

{f*g){x,y) = J2 fi^,^)9{z,y), (2.2) 

z^[x,y] 

where 

[x,y] := {z -.x ^y} (2.3) 
is the segment from x to y. The multiplicative identity is denoted by 5, where 

Six^y) = 1' = (2.4) 
otherwise. 



We have the zeta function 



a^,y) = '^""-^ (2.5) 

' otherwise. 



We think of multiplication by C as 'integration', since 

{C*f)ix,y) = f^^^y)- (2.6) 

The convolution inverse of C, is the Mobius function /i, which satisfies the identity 

Kx,y) = {fJ^*C){x,y) = V n{x,z) = \ ^ (2-7) 

^ , otherwise. 

ze[x,y\ K. 

We will use Mobius inversion on functions from P to C. The incidence algebra acts on 
functions (on the left) as follows. For f E A and (7 : P — t- C, we define 

{f*g){x) = Y.f{x,y)g{y), (2.8) 

and the Mobius inversion formula is given by 

/ = {C*9) ^ 9 = (2.9) 

or, more explicitly, 

(Vx)/(a;) = J^aiy) ^ (Vx)(?(x) = J] i/)/(i/). (2.10) 

y>zx y^x 

The Mobius function of Ii{n) is known (see for example URotl ): if F ^ G and the i-th 
block of G is a union of hi blocks of F, then 

f^{F,G) = {-1Y(^)~^(G) -Q^^^ _ (2.11) 

1=1 

The coefficients //(O, F) and fi{F_, N_) will often show up in our sums and are given by 

HE) 

KQ,F) = (-l)"-(^) - 1)!, (2.12) 

1=1 

f^iZm = (-l)K^)-i(z.(F)-l)!. (2.13) 
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We make extensive use of the following definition. 

Definition 2.1. If F_ :< Ge Il{n) are partitions, we say F_ is a 2-refinement ofG_(or G is 
2-coarser than F) if every block ofG is a union of at most 2 blocks ofF\ If only one block 
decomposes, we say G covers F. 

Covers and 2-refinements arise in our sums, and we note that in these cases the Mobius 
function simplifies to G) = (— Also, if G covers F via the decomposition 



Fi U Fj = Gk, it's easy to see that 








KO,G) 





(2.14) 

More generally, for a 2-refinement F ^ G, let F[ G [F_, G\ be the partition obtained by 
only decomposing the l-th block of G into blocks from F, say Gi = Fi-^ U F/j. Then 

m,F) _ yr KaH) _ ...(F,-.(G)Tj i\Fh\-m\Fi,\-iv. 
Mo,G) \i~m,G) ^ 11 {\Gi\-i)\ ' ^^-'"^ 

where / runs over the blocks Gi that decompose in F. 
Definition 2.2. IfF_ is a 2-refinement ofGj we define the sets 



SiF,G) = {I 

S'{F,G) = {/ 

WiF,G) = {I 

W%F,G) = {I 



Gi decomposes in F}, 

Gi is a block ofF}, 

Fi joins with another block ofF_ in G} , 

Fi remains a block in G}, (2.16) 



soSUS^ = {l,..., u{G)} andWUW = {1,..., z/(F)}. 

Remark 2.3. Given G G Il{n), a 2-refinement F^is uniquely specified by a choice of blocks 
S C {1, . . . , u{G)}, and, for each I E S, a choice of decomposition Gi = Hi U iJf. (If 
|G;| = Ifor some I G S, there are no valid decompositions ofGi.) 

Conversely, given F, a partition G 2-coarser than F is uniquely specified by a choice of 
blocks C {1, . . . , iy{F_)} with \ W\ even, and a way of pairing up the elements ofW. 

2.2. The combinatorial obstruction. We can now clarify some of the obstacles we need 
to address. 

The first reason the random matrix theory and number theory densities in UGaol ap- 
pear different is as follows. In the random matrix theory density, for each partition F = 
{Fi, . . . , Ffc} we at one point consider all ways of decomposing each of the blocks Fj into 
exactly two proper nonempty subsets. That is, we consider all the 2-refinements F^ of F. 
On the number theory side, we instead consider all the ways of pairing up (some or all of) 
the blocks Fj. In other words, we consider all the partitions G of which F is a 2-refinement. 

Because the counting is 'backwards' here, the terms appear very different from those 
encountered on the other side. By reindexing these sums appropriately, we are able to match 
up the parts of the random matrix and number theory densities related to 2-refinements. We 
then reduce the remaining difference to a Fourier transform identity. 
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We verify this remaining Fourier transform identity up to the case n = 7 by breaking 
down the remaining combinatorics. The difference between our approach and Gao's is as 
follows. Gao verified the cases n = 1, 2, 3 by using various ad hoc Fourier Transform 
identities, and explicitly computing formulas for (sums of) integrals over certain regions in 
M" (n < 3), such as (equation 5.11 from LGao D : 

foo pui — l pUi—U2 — l 3 



1il>l+«2+«3 10 

R3,^ n^*("*)^"^ = / / / \[kur)du,, (2.17) 



,_i ^0 ^0 ,_i 



and showed that these sums yielded zero over various sub-regions of the support region 
+ \u2\ + \uj\ < 2. In contrast, we will write 



\\fi{ui)dui = / x{ui-U2-u^)\\fi{u,i)dui, 

i=i -^^io i=i 

/ r3 Ylfi{ui)dui = / x{-ui + U2-U3)Y[fi{ui)dui, (2.18) 



>0 

Ml>l+M2+M3 

3 



where x is the indicator function of the interval [1, oo), and show equality by analyzing the 
combinatorics of various sums of products of indicator functions. 



3. Recasting the expansions 

In this section we rewrite both sides to facilitate the comparison, and reduce the problem 
to a Fourier transform identity. To state the random matrix theory expansion we need the 
following definition. 

Definition 3.1 (x*). For an integer k > 1, the sum of indicator functions xl '■^ 
defined by 

Xl{Ul,...,Uk) = ^ IJJx(M7r(l)H h - M7r(fe)) j , (3.1) 

neSk \i=l J 
7r(l) = l 

where Sk is the symmetric group on {1, . . . ,k}, we sum over the {k — 1)! permutations 
fixing 1, and x is the indicator function of [—1,1]. 

We occasionally refer to xh^ where G is a set (generally a block of a partition G E n(n)); 
the definition is the same as above (with k = \G\) and G is understood the set of indices 
for the variables Ui. 
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3.1. Recasting the random matrix side. The n-level eigenvalue density for USp (see 
equation (4.12) in IIGaoll ) is 



/ n M^w^ri{x)dx = i-^r-"^-^ n ^')' (3.2) 

i=l Gen(n) (=1 



where 



Pi = {\Gi\ [ Giiu)du+ [ Gi{x)dx], (3.3) 
Qi = - J2 - - 1)! / \u\H{u)H''^iu)du, (3.4) 

= ^ / U\Gi\ -1)1- xhii'^h^---^'^i\Gi])jYlfiiui)dui, (3.5) 

with G = {Gi, . . . , G^(G)} and Gi{x) = Hiec, M^)- ^1^°' *e sum Y.[h,h-] ranges over 
the ways of decomposing Gi into two proper, nonempty disjoint subsets H and H'^, and 

-^('") — Ylii^H fi(^) similarly for H''. Except for Lemma [X4l we do not need the 
expansion of "^ritil ^ 

In this section we alter this expression in two ways. First, we rearrange the formula so 
that the Qi terms (involving decompositions of the blocks of G) are put in a form described 
by 2-refinements of G. When we work with the number theory side, we perform a similar 
rearrangement that makes it easy to see the correspondence between these terms. The 
second improvement is to reduce the number of Ri terms we must analyze by showing that 
many (in fact most) of these terms vanish due to support restrictions. 

3.1.1. Reindexing the RMT side. We first recast the above formula in terms of 2-refinements 
of G. 

Lemma 3.1. Equation (13.21 ) is equivalent to 

„ n 2rof 

/ Wmwi^lwx = E2""'^-V(o,z)^(z,G)n(^/ + c^/), (3.6) 

where Y1^f<g ^^^^ over all the 2-refinements ofG (including G_ itself) and 

DiF,G) = n / \u\Hiiu)H[{u)du, (3.7) 



leS{F,G) 

1 

2 



Ai = I Gi{u)du+ I Gi{x)dx, (3.8) 
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and Hi U = Gi is the decomposition ofGi in Fj with Hi{u) = JJ^ fi{u) and similarly 

^ is//, 
for Hf{u) (note empty products are 1 ). 

Proof. We view the sum Y1,[h <I3.2I) as a sum XIf^^g strictly finer partitions 

F_ ~< G_ that are covered by G_ via a decomposition of Gi into if U H^. Note that if Gi 
is a singleton set, then we take the empty sum to be 0. Also, we pull the (IG^I — 1)! and 
(— ) factors to the front, to make a /i(0, F) coefficient. From (12.141 ) we have 

^<i^-0 (3.10) 



(1^1 


-1) 




-1)! 


m 


-1)! 



The new RMT formula is then 



nib " 

/ WfiimtjIWx = 5^ 2'^-^(^V(O,G)n(^^ + ^^ + ^0, (3.11) 



i=l Gen(n) 1=1 

where A; and are as in equations (|3.8I) and (13.91) , and 

where if (m) = JJ^ /^(n), and similarly for ii''=(M). 

Now, we begin expanding the product Yli^i + Ki + Gi) to work directly with the Ki 
term. The goal is to re-express these terms as sums over 2-refinements of G. We have 

i=i Gen(n) ^ sc{i,...,u{G)} les les^ ' 

(3.13) 

We first have the following lemma, which converts the Ki term from a sum over parti- 
tions covered by G_ into a sum over 2-refinements of G_. 

Lemma 3.2. LetG_E Il{n) andletS C {1^ ... ^u{G)} be a fixed subset (i.e., a fixed choice 
of blocks ofG). Then 

2ref,5 

MO,G)n^/ = E/^(^'^)^(^'^)' (3.14) 

leS F-<G 

where Ylip'^G ^^^^ ^^^^ 2-refinements F of G_ such that S{F^,G) = S is the set of 

blocks ofG that decompose in F_. The term D{F^, G_) is as in (13.71) and Ki is as in (|3.12l) . 

Remark 3.3. In order to have any 2-refinements F^of G_ in the right-hand side of (13.141) 
above, each of the blocks Gi (I G S) must not be a singleton set. Note (13.141) holds either 
way. If Gi is a singleton set for some I G 5", the Ki factor on the left-hand side and the 
entire right-hand side are both empty sums, hence zero. 
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Proof of Lemma \3?2\ Expanding the left-hand side, we have 

m.G)\{Ki = MO,G)n E / HHiu)H^{u)du. (3.15) 

les les F^G — '—^ 

When we expand this sum, we obtain a sum of terms, each of the form 

m,G) ■ n^^^^ / HHi{u)H!{u)du, (3.16) 

where is the partition covered by G by decomposing the block Gi = HiUH[ and leaving 
the other blocks of G_ unchanged. 

Let F ^ G be the partition obtained by decomposing all the Gi this way. Then each 
summand corresponds to a unique such F, a 2-refinement of G with S{F_, G) = S. By the 
identity (|2.15l) . the fi coefficient becomes 

''<^'^> ■ - ^(aC) . ^ ^ ,(0,£), (3.17) 

so the term simplifies to 

/^(0,£:)n / \u\H^iiu)H[iu)du = f^{0,F)D{F,G), (3.18) 

as desired. 

Conversely, every 2-refinement G with S{F_, G) = S arises (once) this way, so the 
two sides of (|3.14l) match. □ 

We now return to the proof of Lemma 13. 1[ Next, when we sum (|3.14l) over all S C 
{1, . . . , i^{G)}, we get a sum over all the 2-refinements F of G (including G itself, from 
5 = 0). We have 

„ n 2ref 

i=l G£U(n) F^G «e5<= 



GeU{n) 

2-refine] 

of Lemma im □ 



where YI'f'^g ^^^^ 2-refinements of G (including G itself), completing the proof 



3.1.2. Expanding the Gi terms. We expand and simplify the H/l^^ + ^i) term. The fol- 
lowing lemma drastically reduces the number of terms we have to analyze. 

Lemma 3.4. Let Gi and Gk be disjoint subsets o/{l, . . . , n}. Then 

^' ■ ^ L - al&) n ■ L - (Ift^) n S'^^^- 

= 0, (3.20) 
where xhi (^) '■^ shorthand for Xgi (^n ; • • • ; '^i\Gi\ )' defined in equation (13.11) . 
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Proof. Since Gi and Gk are disjoint, we must have either 

5^supp(/i)<l or 5^ supp(/,) < 1, (3.21) 

since the total support is less than 2. Without loss of generality, assume G^'s total support 
is less than 1. Then 

le^nMiiH \-ei^Ui^ \ < y^\ui\ < 1 (3.22) 

^ ■' ^ — ^ 

in the region of support, so xC12fi ~ ^ ^^^^ ^' ~ ^^^^^ ^hi ^ °f 
{\Gi\ — 1)! products of x's, the Gi integrand is identically 0. □ 

To emphasize the significance of this lemma, we note that instead of having to expand a 
product of the form nf=i(^« + ^i) i^^to 2^ terms 

k 

YliA + Ci) = n^'H^'' (3.23) 

1=1 UC{l,...,n} IGU l(f.U 

we only end up with k + \ nonvanishing terms: 

k k k 

\[{Ai + Q) = n + E • n (3.24) 

1=1 1=1 1=1 i'=ii 

Combining Lemmas [3 .41 and [3711 yields the following. 
Lemma 3.5. With notation as in Lemma 13.71 (|3.2I) is equivalent to 



\{mwi^i{x)dx = 

i=l 

2rcf . . 

E E2"~'^-^Mo,F)D(F,G) n^'+E^'-n^'' • ^^-^^^ 

Gen(n)ZdG ^IGS'' I'^l ' 

The expression (13.251) is the one we use when we start matching terms with the number 
theory (NT) side. 



3.2. Recasting the NT formula. We now recast the NT density as a sum over 2-refinements 
of partitions, bringing it closer to the RMT formula established in Lemma [XT] This allows 
us to fully match one set of terms appearing on both sides. We then alter each formula 
slightly to reduce the problem to a Fourier transform identity, relating the terms Ci on the 
RMT side (equation (13.91 )) to the integrals over ]R>q on the number theory side (equation 
(I33T]) ). 

Gao's expression for the n-level density of zeros of quadratic Dirichlet L-functions, 
which we abbreviate as IVg"-*, is (adapted from equation (2.16) in UGaol ): 
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2 '^(^) 



llmW^''\x)dx = Jim — E 2--'^^^^fi{0,F)l[{A + Bi), 

i=l deD{X)Fen{n) 1=1 

(3.26) 



Ai = I Fi{x)dx - ]- / Fi{u)du, (3.27) 



where 



2 

Here o? is the conductor, F_ = {Fi, . . . , F,^(^f)} and F;(x) = YlieFi Zip o^^^ the 

primes and (y) is the Legendre symbol. 

Note that the Ai terms are independent of d and X. Hence, if we expand the products, the 
Ai terms can be pulled past limx^oo J2deD{x)^ making their contributions easy to analyze: 

deD{X) ^ I ^ ^ V ' ^ I ' ^ d£D{X) v ' 

The main difficulty comes from the expressions 

lim V WBi, (3.30) 

d£D{X) lew 

where C {1, . . . , i^(-F)}, since the Legendre symbol (y) in the series (|3.28l) introduces 
a dependence on d and X. 
For these, Gao develops the following formula (see equation (3.13) in nOaoH ): 

Lemma 3.6. Let F = {Fi, . . . , -Fi/(£)} be as above, and let W O {I, . . . , I'iF)}. Then 
lim — V Y\ Bi = (3.31) 

d€D(X) lew 

n + A^iH^i^ jj r ..F.,(..)F,(..)rf.. + (-2)1^1- E (E n 

^ ^ {A;B) i=l W2^W ^ (C,D) i=l 

jVK2| even 

M,F,^(n,)F,^(n,)rfM,') • i f^^^ ( J] (_l)l^l^( J] _ J] u,)] H^.l^O^^^), 

/ V "'K>(,2 V IGW^ 7= I ' W^ ^ 

where W2 = W \ W2, and the notations J2{a-b) Y^(cd) ^^^^ ways of pairing 
up the elements ofW and W2, respectively. Also, x is the indicator function of the interval 
(1, 00). Empty products are 1. 
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We now obtain the general formula for the n-level density by combining expressions 
(13.261 ) and (13.311 ) and using the expansion 

HE) 

i=l WC{l,...,iy{F)} W 

3.2.1. Reindexing the NT side. We put the NT formula in a form closer to the RMT for- 
mula, as a sum indexed by 2-refinements of partitions. We establish the following. 

Lemma 3.7. Gao's expression (|3. 261) /or the NT density is equivalent to 

n 

l[Mx)W^'\x)dx = 

i=l 

2ref . . 

j2 E r-'^'^^^m,F)D{F,G) n - 2 E ^(^' ^) n ^0 ' ^^-^^^ 

where ^Y^p^^Q runs over the 2-refinements F ofG_, the sets S{F_, G) and S'^ = S'^{F, G) are 
as in Definition \2.2\ and 

D{F,G) = n / \u\Mu)H[iu)du, (3.34) 



leS{F,G) ■ 
-1 

y 



Ai = ^ I Gi{u)du+ I Gi{x)dx, (3.35) 



E{G,T) = 2l^l f^^ ('^(-i)l^l^(^u. - E«^)) U^Mdui, (3.36) 

where for I G S{F_, G), Gi = Hi U Hf is the decomposition of the block Gi into blocks 
of Fj and x is the indicator function o/(l, oo). Empty products are 1 and empty sums, in 
particular E{G, 0), are 0. 

In order to prove Lemma [3771 we first alter formula (13.311) for the Y[ Bi terms. 

Lemma 3.8. Let F_ E n(n) and W <^ {1, . . . , i^{F_)}. The following formula is equivalent 
to Gao 's Lemma \3.31[ 

x'-^oo AX 



d£D{X) lew ^ ^ ^ GtF ^ 

, 2cor,VK2 X 

2 E ( E 2'^^-^''^^-^^(£.G))E(F,iy2'), (3.37) 



2 

W2CW ^ gyp 

\W2 \ even 

where X]g°f^ X]g°f^^ ^^^^ partitions G 2-coarser than F_ with W{F_, G) 
W and W2, respectively, and the other notation is as in Lemma \3.7\ 
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\W\I2 2cor,VK 

2"^'En/ y'KXn)\{u)du = ^ 2^(^)-'^(^)D(F,G), (3.38) 



Proof. If |iy I is even, we claim 

\W\I2 2cor,VK 

where the left-hand side is as in Gao's formula, (I3.31I ). (The same identity holds with W 
replaced by 1^2-) To see this, first note that by Remark [2l3l each way of pairing up the 
elements of an even subset 11/ C {1, ... , corresponds to a unique partition G_ that is 

2-coarser than F, with G) = W. This correspondence is one of the key ingredients, 

as it allows us to begin expressing the sum in terms of 2-coarser partitions. Later we will 
switch orders of summation, converting sums over 2-coarser partitions into sums over 2- 
refinements, which is what we have on the RMT side. 
Thus 



\W\/2 2cor,W\W\/2 

2"^' Z n / uUuAiu)du = 2l'^l ^ H uFaAu)hAu)du. (3.39) 



For the integrands, observe that pairing Fi^ with Fi^ to form a block Gi of G_ is equivalent 
to decomposing Gi into two subsets Gi = Hi U Hf with Hi = Fi^ and Hf = Fi^ . Since 
each Fi- is an even function, we can replace with ^ J^, and u Fi^{u)Fi^{u) becomes 
\u\Fi^ {u)Fi2 (m) , as in the definition of the term D (F, G) : 

2cor,Ty \W\/2 2cor,Vy 

2'^' E n / uFMFMdu = 2l^l 5^ 2-l^l/2D(F,G). (3.40) 

GtE i=l G>F 

Finally, for the 2'^'/^ coefficient, observe that ^{F) — ^{G) = \W\/2, since each pairing 
reduces the total number of blocks by 1. We apply identity (13.381) to both the 'Y1i{a-b) 
X](C D) terms in Gao's expression (13.311) to obtain the desired form (13.371) . □ 

We return to the proof of Lemma [3771 Applying Lemma [X8] to Gao's expression (13.261 ) for 
the NT density gives 

n 

\[h{x)wi-\x)dx 

i=l 

Fen{n) VKc{i,...,y(F)} i<$_w deD{X) lew 

Si + S2, (3.41) 
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where 



(2cor,W ^ 

|W| even 

(3.42) 



_1 2»-e)M0,F) 

Fen(n) WC{l,...,i/(F)} 10V 

(2cor,W2 X 
^ 2-^^^-'^^^D{F_,G)\e{F_,W^). (3.43) 
n-^TT / 



IW2I even 



We work with and S2 separately, since S2 includes an extra summation. For 5"!, we 
have 



, 2cor,VF ^ 

Fen(n) WQ{l,...,u{F)} l^W ^ GhZ ^ 

\W\ even 

2cor,W 

= E E E 2""''^-^/^(0,F)D(F,G)n^/- (3.44) 

Fen(n) VKC{1,...,!/(F)} GbZ l^W 
\W\ even 



2cor,PF 

Now the double sum is equivalent to summing over all the partitions 

W<Z{\,...,v{F)} QhF 
\W\ even 

G_ that are 2-coarser than F_, i.e., X^^^'p, since every such G_ arises exactly once this way 
(including G = F, from the case W = 0). Note that the set {/ ^ W} is just the list of 
blocks that are common to both F_ and G, so it is the same as S'^{F, G). By switching the 
order of summation of F_ and G, we obtain 



2cor 

^1 = E E2""'^^-V(o,£:)/^(z,G) n 

2rcf 

= E E2""''^-^)"(^'^)^(£'^) n (3.45) 

Gen(n) ZdG leS'={F,G) 
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The term 5*2 is more delicate, since it includes an extra summation: 



Fen{n) WC{l,...,u(F)} IfW 

(2cor,W2 \ 
n^TP / 



W2CW ^ GtF 
IW2I even 

2cotW2 

= -2 E E EE 2-^(^^,,{0,F)DiF,G)E{F,W\W,)llAi. 

Fen(n) WC{l,...,u(F_)} W2CW QtF IfW 
\W2 \ even 

(3.46) 

We switch the choice of subsets W, W2 Q {1, . . . , i^(F)}- In particular, the choices of 
W and W2 effectively partition {1, . . . , i^{F)} into three disjoint subsets: 



W2 (with I W2 1 even) 
W-W2 
{l,...,uiF)}-W 



lists the blocks of F merged to form blocks of G in the D{F, G) term 
lists the blocks to go in the E term, 
lists the blocks to go in the Y[ -^i term. 

We switch this so that W2 is chosen first, which allows us to pull the X]f<g^^ 
front. In other words, we choose W2, followed by a disjoint set T C {1, . . . , v{F)} — W2, 
which lists the blocks to go in the E term. With this change, the E{F_, W \ W2) is replaced 
by E{F_, T), and Hi^w -^i becomes n/^vF2UT -^i- P*^^^ the Y1^g>^f^^ outward: 



2cor,VV'' 



^2 = E E EE 2— (^)/.(o,F)i^(F,G)£;(F,r) n 

£en(n) W2C{1,. ..,;/(£)} rCVF| G^F l^TUW2 
IW2I even 

2cor,Vy2 

= -2 E E EE '2''-'^^^m,F)D{F,G)E{F,T) J] ^^ 

Fen{n)W2C{l,...,u{F)} GtF TCW^ lfTUW2 
IW2I even 

(3.47) 

Now the summation X]|iy2| even is the same as what we encountered in our analysis 

of iSi, a sum over all partitions G that are 2-coarser than F_ (including G_ = F_, from the 
case W2 = 0). The set W2 is the same as 5''^(F, G), the list of blocks common to both 
partitions, so we rewrite ^j^cv^i Stc5'=(f g) ^^'^ Y[i^tuW2 ^^^^ Yiies^-T- obtain 



^ 2cor 

= -2 E E2""'^-^M0,Z)i^(Z,G) Yl n (3.48) 



£en(n) GtF TCS''{F,G) l&S^-T 
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Now we switch the and converting S2 into a sum over 2-refinements: 



2rcf 



Gen(n)F^G TC5=(F,G) leS^-T 

Finally, we rewrite E{F,T) = E{G,T). This is just a relabeling, since T C 5*^ the 
set of blocks Fi e F that are unchanged in G_, and the integral over M>o in the definition 

of E (equation (|3.36l) ) only involves the functions Fi{ui) where / G T. Nonetheless, it is 
important as it expresses the E term in terms of the outermost summation 

Putting together our expressions (13.451) for Si and (13.491 ) for S2 yields an NT formula 
expressed in terms of 2-refinements: 



\{h{x)Wl^\x)dx = 

i=l 

2ref . . 

This completes the proof of Lemma [X71 □ 

Remark 3.9. The key step in the proof of Lemma WJ\ was to switch the orders of summation 
of F_ and G in the NT density, replacing a sum over 2-coarser partitions by a sum over 2- 
refinements. On the RMT side, 2-refinements already appeared naturally as products of the 
terms Qi (equation (|3.4I) ). and thus no switch was necessary. 

3.3. Reducing to the Fourier identity. Lemmas [3.5l and [X7] establish the following forms 

(n) (n) 

for the RMT and NT density expressions, VFygp and Wq . 



RMT: / Y\f^{x)wlj%{x)dx = 

2ref . . 

E E2"^''^-^M0,F)D(F,G)( n^ + l^a-n^'' ' (3.51) 

Gmin)F^G ^«G5= leS'' I'^l ' 

„ n 

NT: / T\f,{x)W'^^\x)dx = 

2rcf , . 
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where J2f<g ^^^^ ^he 2-refinements F of G, and 

D{F,G) = n / HMu)H[{u)du, (3.53) 

1<^S{F,G) '^^ 

Ai = ^ I G'i{u)du+ I Gi{x)dx, (3.54) 
Q = -[ fl - TT kn.)du,, (3.55) 

JR-^a \ ICT I / 

where for / G S{F^,G), Gi = HiU Hf is the decomposition of the block Gi into blocks 
of F, and x is the indicator function of (1, oo). Empty products are 1 and empty sums, in 
particular E{G_^ 0), are 0. 

We have some cancelation right away: namely, the terms n«G5'= without Gi or E{G_, T) 
factors match, since the new expressions count all the D{F^,G) factors the same way on 
both sides. Compare this with the original density expressions (|3.2I) and (13.261) . which only 
make it easy to see equality between the terms with all Ai factors (without any of the C/, 
E{G,T) or D{F_,G) factors). Those terms show up in the new expressions as the trivial 
2-refinements where F_ = G_. 

Unfortunately, with the expressions above, the sums do not match term-by-term: the E 
terms combine across many different 2-refinement pairs (G, F). The goal of this section is 
to reduce the Density Conjecture to an identity relating the E{G,T) to the Gi terms. We 
use Mobius inversion to express the identity in a fairly simple way. We then verify the 
identity for n < 7 by breaking down the remaining combinatorics. 

3.3.1. Isolating the Gi and E terms. By canceling the matching terms in the two densities, 
we are reduced to showing equality between 



2ref 



RMT = 5^ 5^2--(^^)/i(0,F)D(F,G)EQ-n^z', (3-57) 

Gen(n) F<G «G5= I'^l 



2ref 



NT = -- E 5^2"-'^(SV(0,F)2;(F,G) 5^F(G,T) J] (3-58) 

G€U{n)F<G TCS'= IdS^-T 

with notation as in (I3.53l) - (|3.56l) . Note that these expression are not the same as the n-level 
density expressions, (|3.51l) and (13.521) : all the matching terms have been removed. 

In this section, we rewrite the G and E terms to depend only on T, not G_. This allows 
us to pull them outside the summation We show the following. 
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Lemma 3.10. Equations (13.571 ) and (13.581) are equivalent to 
RMT = \ Yl (2l"IC(0^)) ■RestlW'^), 

WC|l,...,n| ^ ^ 



WC{l,...,n} ^ T&n{u) 



(3.59) 
(3.60) 



Rest(W") 
C{T) 

E{T) 



2ref 



^ ^2l"^l-(^V(0^^,F)D(F,G)nA 



l&S" 



(3.61) 



(3.62) 



E (-l)'"+'x(E"«-E"«)) n'^("')^«'' (3.63) 



an J '■^ minimal element ofIl{U) (all singleton blocks), and the rest of the nota- 
tion is as above. Note that fi{T,N) = - 1)' ■^o fi{0^,Ku) = 

{-ir\-\\u\~i)\. 



Proof. We first work with the NT side. First of all, from the definition in (l336l) . E{G, T) 
is an integral involving only the functions 



G,{u) 



u) 



(3.64) 



from the blocks Gu I G T. We obtain the E term by choosing a partition G G n(n), 
followed by a 2-refinement F, followed by a choice of blocks T C S''^(F, G) . 

To isolate the E term, we switch orders. We first choose a subset W C {1, ... n} of test 
functions and a partition T G n(W), and then choose a partition of the remaining elements, 
G' G n(W"), and a 2-refinement F' of G'. We use E{T), as defined in (13:63]) . Note that 
E{G,T) = (-1) ■ 2'^(^)F(7:). (Compare (13:63]) and (1X56]) .) 

So, for each term we have 



E{G,T)D{F,G) H A = -2''^^'>E{r)D{F:,G!)l[Ai 



(3.65) 



/gS'=(f,g)-t 
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For the 2"-'^(^V(0, Z) coefficient, we liave to pull out a factor of 2l"l-'^(^V(£2w,2I) 
(where is the minimal element oiYi{U)): 

HE) 

2"-'^(^)/i(0,F) = 2"-'^(^)(-l)"-"(^) - 1)! 

1=1 

= (2l"l-^(S^(o^,r))(2l"^l-(^V(0^^,Z')). (3.66) 

Note that the 2"'"® will cancel with the 2^"® coefficient on £"(71) in (13.651) . 
This gives the desired expression for the number theory side: 

UC{l,...,n} ^T£U{U) ' G&l{W)F_-<G leS^ 

(3.67) 

as desired. 

We proceed similarly for the random matrix theory side. First of all, the Ci term ap- 
pearing in (|3.55l) always has the test functions arranged as HiGG; fii'^i)dui, for some block 
Gi of G, with each test function Fourier-transformed separately. Thus with notation as in 
(13.621) . the Ci term always takes the form 

where U = Gi C {1, . . . , n} and is the minimal element of n(W), with each test 
function in its own (singleton) block. (It is nonetheless necessary to define C(T) for any 
partition T G n(W), in order to use Mobius inversion later.) 

The argument is now similar to (in fact more straightforward than) the NT side. The Gi 
term in (13.581) arises choosing a partition G E n(n), a 2-refinement F_ of G, and a single 
block Gi G S^IF, G) to put in the Gi term. 

We switch orders. We first choose a subset U C {1, . . . ,n} (from which to get a G{0^) 
term), then choose a partition G' G Il(U^) and a 2-refinement F' of G'. As with the 
RMT side, the D(F^,G) and Ai terms are unaffected, but we have to break up the Mobius 
coefficient n{0, F_). As the — 1)! factor cancels the ^i^i-*^,^^, factor on the G{Ojj) term, 
the coefficient becomes 

(^\\U\-l 

2— (^V(0,F)-Q = 2— (g)(-l)— (^) n(|F,| - 1)! ■ ^^f {Q^) 

= (2l"l-^G(0^)) (2l"1-(^V(0^.,Z')) • (3.69) 
So for a single term, we have 

2--'^'>fi{0,F)D{F,G)Gi l[A, = 2l"l-iG(0^) [2l"^l-'^(^V(0i,c, F')/^(F', G') J] 

(3.70) 
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Thus the RMT side becomes 

^ 2rcf 

RMT = - Yl 2l"IC(0^) E E2'"''^'^-V(0,F)D(F,G)nA, (3.71) 

UC{l,...,n} Gen(W=)F^G «G5<: 

which is the desired expression. □ 

3.3.2. The Fourier identity. Lemma [3.101 reduces the density conjecture to showing that 
the two expressions (13.591) and (13.601) are equal. These expressions separate the C and E 
terms from the others, so the question is: how do they match up? We believe the follow- 
ing conjecture, which essentially says that they match term-by-term in the form given by 
Lemma [3. 101 

Conjecture 3.11 (Fourier Identity 1). With notation as in Lemma \3.10\ and * denoting 
(incidence algebra) convolution, C is the Mobius transform of E: 

C = fi* E, or equivalently C = E, (3.72) 

as functions onIl(U). 

In particular, the identity we need, which we apply once for each subset U C {1, . . . , n}, 
is simply 

Conjecture 3.12 (Fourier Identity 2). With notation as in Lemma 13.701 

C{0^) = J2 l'iQ^^I-)E{r). (3.73) 
ren(w) 

Equivalently, by Mobius inversion, 

CiT) = E{0^). (3.74) 

Teu{u) 

It is clear that Conjecture [3TTT] implies Conjecture 13.121 and Conjecture 13.121 implies the 
Density Conjecture. In fact, Conjecture l3.12l is equivalent to Conjecture l3.1 IK and. as such, 
it must be true by the results of HER-GRL though there should be a purely combinatorial 
proof of this fact). The equivalence stems from the fact that if F G n(?7,) has k blocks 
Fi,...,Fk, then C{F) and E{F) are the same as the integrals C(0), E{0) for O e U{k), 
using the Fi(x) as a new set of test functions. Thus, if the identity (|3.74l) holds for k < n 
(for all choices of test function), so does the identity (13.731) . 

For our purposes, identity (13.741 ) is preferable since, in contrast to (13.731) , all the sum- 
mands are easy to convert to integrals over the same region (]R>o), which we do in ^ We 
summarize our results so far. 

Theorem 3.13 (Reduction to Fourier Identity). Let fi, . . . , fn be even test functions with 
fi, . . . , fn supported in J2^=i l^il < 2. The Fourier identity (13.741 ). with U = {1, . . . , n}, 
implies the Density Conjecture for quadratic Dirichlet L-functions, 

/ Uf^^^'^^uspi^^dx = / llf^{x)wi^''\x)dx. (3.75) 
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In the next section, we study the Fourier identity further and put it in a 'canonical' form, 
which we use to verify the cases n < 7. We prove 

Theorem 3.14 (Density Conjecture, n < 7). With notation and assumptions as in Theorem 
\3.13\ the Fourier Identity (13.741) holds for n < 7 (for allU C {1, . . . , n}j. 

The above immediately implies our main result. Theorem II. 1[ In particular, for < 7 
the n-level density of zeros of quadratic Dirichlet L-functions {L{s, Xsd)} (with d G N odd 
and square-free) is the same as the n-level eigenvalue density of the Unitary Symplectic 
Ensemble (for test functions where the sum of the supports is at most 2). 

The remainder of the paper is devoted to the proof of the Fourier Identity and Theorem 

EH 



4. The Fourier Identity 
In this section we consider the Fourier Identity (13.741) for a fixed n: 

Fen(n) 

where for F G n(n), 

1 r f V^^^ ^ 

C{F) = - / MZ,iV) + (-l)^(=^)x:(F)(wi,---,MZ)) l\Fi{ui)dui, (4.2) 

^K5o V/C{l,...,n} / ^ 2=1 

see Definition 13. II for the definition of xi^(F)- 

We first reduce the Fourier Identity to a canonical form. Our method is to convert all 
the summands C(F) to integrals over ]R"q and reduce to a sum of products of indicator 
functions of the form 

^(eiMi H VSnUn), X = I(l,oo), (4.4) 

with each Si = ±1. We examine, for each term, the set : = +1}, and use combinato- 
rial arguments and the assumption ^ supp(/j) < 2 to simplify some terms and show that 
others are identically zero. We start by reducing to a canonical form: 

Proposition 4.1 (Fourier identity, canonical form). With notation as above, 

„ n 

E C{F) = / Co(wi,...,M„) J]7,(M,)rfM„ (4.5) 

F£U{n) ''^h i=l 

where 

k 

Co{u,,...,Un) = E E i-^r''U(XAAJ + XA^Aj), (4.6) 

JC{2,...,nj chains A i=l 
'Ai = {1} 
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where A ranges over the chains Ai d ■ ■ ■ C A^, such that Ai = {1}, and for a subset 
C {1, . . . , n}, we write 

Xw-= x{^Ui-^Ui), X = I(i,oo), (4.7) 

and A denotes symmetric difference. 

Hence, comparing (|4.3I) and (14.51) . the identity of functions 

CoK,...,n„) = (4.8) 

/C{l,...,n} 

on the simplex \ui > 0; ^ Uj < 2} C M" implies the Fourier identity (13.741) . 

While we are unable to prove (14.81) for all n, we give a method for checking it for specific 
values of n. Our method is partly ad hoc, but suffices for n < 7 (for larger n, the computa- 



tions are the same, but become intractable). We note that for n = 1, 2, 3, the identities (|4.1I) 
appear in UGaoH as, respectively, the last (unnumbered) equations on pages 57 and 58, and 
equation (5.7). 



4. 1 . The canonical form. We first reduce (14.11) to a sum of indicator functions by convert- 
ing all integrals to the region ]R>o. 

Lemma 4.2. With notation as above, the Fourier Identity (14.11) for n> \ follows from the 
equality of indicator functions 

J2 C\F) = E\0) (4.9) 
£en{n) 

on the region {0 < < supp(/j), i = 1, . . . , n)} C ]R>o, where 



(4.10) 



^'(^) = \Y. (/^(z,iv) + (-i)'^(^)x:(^)(5^e.«.,..., ^^«^)): 

i=l,...,n 

E\0) = Yl (4.11) 

/C{l,...,n} / I" 

(Note that E'{0) is just the sum of indicator functions in the integrand ofE{0_).) 

Proof. Our test functions /i, . . . , /„ are all even, so we have two identities. First, for any 
partition F G Il{n) having k blocks, 

/ g{ui,...,Uk)Y\^i{ui)dui= / 5((V'ui,..., Ui)Y\fi{ui)dui (4.12) 

■^K* t-^ >^R" ,-cP 

holds for any integrand g via a linear change of variables. Second, 

/ . . . ,M„) ]T/i(Mi)(iMi = / . . . ,£„M„) ]T/i(Mi)(iMi (4.13) 

i=l .,=±1 -^K>o i=l 



i=l,...,n 
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holds for any h. Applying these transformations gives, with notation as in (14. 10k 

„ n 

C{F) = / C\F)-T\Mu,)du„ (4.14) 

■^^>o i=l 

which is the desired expression. □ 
We use the identity 

Xiu) = l-xiu)-xi-u) (4.15) 

to rewrite the left-hand side as a sum of products of terms of the form x{Yl ^i'^i) with each 
Ei = ±1. The advantage of using x throughout comes from not being an even function: to 
know whether xC^ ^iUi) = 0, we need to consider both ^ EiUi > 1 and ^ EiUi < — 1, but 
with X only the first case matters. This will facilitate several simplifications. 

Definition 4.3 (Combinatorial notation). We adopt the following notation: given a term 
x{siUi + ■ • ■ + EnUn) , let A be the set of indices for which Ei = +1 and the set for which 
Ei = —\. We define 

XA ■■= Xi^Ua, - ^UaJ. (4.16) 

A A<= 

This notation reduces arguments about products of x{J2i to combinatorial argu- 
ments about subsets AC {1, . . . , n}. The x* integrand for a partition F_ e n(n) is thus, in 
combinatorial notation, 

7r(l) = l 

where S^(^f) is the group of permutations of {1, ... , z^(F)}. 

Changing the signs of some of the Ei in xa is equivalent to taking a symmetric difference, 
replacing xa with Xaaj, where J C {1, . . . , n} is the set of indices whose signs have been 
changed. We note that yu(F,iV) = (— l)*^^— ^~^(z^(F) — 1)! is the same as the number of 
permutations on the inner sum. We thus write 



C'{F) = 1 J2 Uf,N) 

JC{l,...,n} ^ 

^(^_iy{F) ^ JJ (l-X(F,(i)U-UF,(,))AJ-X(F,(l)U-UF,(i))AJ-) j 

7rGSj,(_F) i=l 
7r(l) = l 

5 E E 



2 

JC{l,...,n} ttGS^^p^ 
7r(l) = l 

HE) 



+ n (X(F,(i)U-UF,„)AJ + X(F,(i)U-UF,(,))AJ- - l) ) • (4.18) 
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In combinatorial notation, the right-hand side of (14.91) is just 

E'm = Yl (4.19) 

/C{l,...,n} 

To shorten the notation, we write the summands in terms of chains rather than partitions. 
Given a partition F_ G n(?7,) and a permutation tt e S,y(^F), such that 7r(l) = 1, we obtain an 
ascending chain 

i^.(i) C U F,(2) C • ■ ■ C U • ■ ■ U F^^k) = {1, • • • , n}. (4.20) 

Thus each choice of F_ and tt corresponds uniquely to a strictly ascending chain 

A: A^CA^C-.-CAk (4.21) 

of subsets of {1, ... , n} such that 1 E Ai and A^ = {I, . . . , n}. The corresponding product 
of indicator functions is then 

(XAiAJ + XAIAJ - 1){XA2AJ + XA-AJ - 1) ■ ■ ■ {XAkAJ + Xa-aj " !)• (4.22) 

Thus we can write the left-hand side of (14. 9|) as 

E ^'(^) = ^ E E ((-l)'"' + n(XA,Aj + X^.Ajc-l)Y (4.23) 

£Gn(n) JC{l,...,n} chains yt,leAi ^ i=l ^ 

where A ranges over all ascending chains of subsets of {1, ... , n} such that 1 E Ai and 
Ak = {I, ■ ■ ■ ,n} is the last (largest) set in the chain. 

Observe that a given product of x terms occurs many times in the sum (14.231 ) when the 
(—1) factors are expanded. We account for this cancelation below and give the canonical 
form of the Fourier Identity. We employ the following standard facts about chains. 

Proposition 4.4 (Sums over chains). Let A and B be chains. If A = Ai C ■ ■ ■ C Ak, write 
k = \A\, and ifB is a subchain of A, write A^ B. For a fixed B with 1 E Bi, 



J2 (-1)1-1 = ^(-1)" = (4.24) 

AtmeAi. I ° otherwise. 



A-={l,-,n} 

See appendix |B] for a proof. 
Lemma 4.5 (Fourier identity, canonical form). With notation as above, 

E cm = E E (-1)""'-" 

Fen(n) JC|2....,n| chains .4 i=l 



{Xa.aj + Xa-aj), (4.25) 



where A ranges over the chains Ai C ■ ■ ■ C A^ (k = \A\) such that Ai = {1}. (We do not 
require = {1, . . . , n}.) 
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Proof. Consider the expansion for ^^gn(n) ^'(E) in <I4.23I ). Since the right hand side is 
invariant under interchanging J and J'^, we replace | X]jc{i n} "^^^^ jc{2 n}' ^^^^ 
we may assume without loss of generality that 1 ^ J. 

Now we expand the (—1) factors in (14.231) . The (—1)'^ cancels with the (—1)''^^, so we 
are left with 

£en{n) JC{2,....n} chains ^,leAi Vl/C{l,...,t} i&W 

Afe = {l,...,n} l^^0 

= E E E (-l)'^'-'"'nfe.A. + XB.A.O, (4.26) 

JC{2,...,n} chains ^,leyli chains B^y4 B 
Ak = {l,...,n} 

where B ranges over the subchains of A (excluding the 'empty chain' with no sets). We 
switch orders of summation on B and A. We have 



chaias A,leAi 0^B^A B AhB,l&Ai, 

Ak = {l,...,n} leBi Afe = {l,...,n} 



and so 



E^'(£)= E E (-i)'"nfe^A.+XB.A.o( E (-!)'■ 

Fenfn) JC{2....,n| chains B B ^>-i3,leAi, 

= E E (-l)""""nfe«^-^ + *^»^^^) (4.28) 

JC{2,...,n| chains B B 

by Proposition |4]4] □ 

4.2. Breaking down the combinatorics. We describe our approach to confirm the Fourier 
identity (14.91 ) in the cases n < 7. These arguments are impractical to do by hand for ri > 4; 
we ran them in Mathematica with code available at 



htt p : / / www-personal . umich . edu/~ jakelev/ 



• http : / /web . w illiams . edu/Mat hemat i cs/s jmiller/public_html /math/papers/ jakel/Fourierldentity . tar| 

The simplifications we use are as follows. 

Lemma 4.6 (Simplifications). Let A, B C {1, . . . , n}. Then 

Xa- Xb = Xa whenever A C B. (4.29) 

Let Ai, . . . ,Ak C {1, . . . , n}, with k > 2. For each i G {1, . . . , n}, let Cj be the number of 
the Aj 's that contain i. Then 

XAi • • • XAfc = if < |A; for each i e {I, . . . ,n}. (4.30) 

The first equation says, equivalently, that given a product xai ■ ■ ■ Xa^ , we need only keep 
the 's for which the subsets Aj C {1, . . . , n}'s are minimal with respect to containment, 
i.e., the x terms having few positive signs. The identity is essentially a formal sum of 
antichains, with additional relations such as (14.301 ). We remark also that for k = 2, ( 14.301) 
is just the statement Xa ■ Xb = ^ when An B = 0. 



n-LEVEL DENSITffiS OF DIRICHLET L-FUNCTIONS 27 

Proof of (14.291 ). li A C B, then B'^ C A'^, so we have the inequalities 

^Ub^ > ^Ua^ and ^Ua^ > ("^-31) 

B A A'' B'' 

Combining these yields that whenever — Y1,a-= > 1' we also have 

^^h-^'^h > ^Ua,-^Ua^ > 1- (4.32) 

So if xa = 1, it follows that Xb = I (if Xa = 0, then both sides of (14.291) are 0). □ 

Proof of (14.301 ). Add the inequalities XlaeA^a ~ Z]a'GA= ^a' > -'- together. Then +Ui 
occurs Cj times and —Ui occurs k — ci times, so the result is 

{2ei - k)ui + ■ ■ ■ + (2e„ - k)un > k. (4.33) 

The condition Cj < |/c is the same as | > 2ej — k, yielding 

+ ■■■ + «„) > k, (4.34) 

that is, Ml + h M„ > 2, violating the support restriction. □ 

Remark 4.7. For k < A, the second condition (14.301) is equivalent to H^Li = ^- F'or 
k > Ait is a stronger condition. 



In sum, our method of verifying the Fourier identity (14.91 ) is to apply the simplifications 
above to the sum (14.251) to simplify and remove terms. We are able to verify the cases n < 7 
this way; forn = 8 the verification becomes intractable, since the number of terms in the 
left-hand side of (14.91) becomes large (2"- sequence A027882 in the Online Encylopedia of 
Integer Sequences). 

Although we cannot prove the identity for all n, we give one last conjecture that indicates 
one way of grouping terms in the identity. 

Conjecture 4.8. For fixed n and J C {2, . . . , n}, let 

^AC{J} ^ ^B:JCBC{2,...,n} ^ 

simp^H = (-1)'^'XAU{1}. (4.35) 

AC{2,...,n} 

Then the inner sum of the Fourier identity (|4.25l) is 

( A\^-\^\Y\(^ N _ /simpi(J;n) J ^ {2, . . . , n} 

>, \-^) I {Xa^aj+Xa-aj) - < . ( T \ , ■ ( \ J fn \ 



chains A i=l 
Ai = {1} 



(4.36) 



It is easy to see, by summing over J C {2, . . . ,n}, that this conjecture implies the 
Fourier identity. The identity is easily checked for J = 0, {2}, {2, . . . , n}; for the remain- 
ing cases, it is sufficient (by relabeling) to consider J = {2, . . . ,i} for 3 < i < n — 1, but 
we do not as yet have a proof. 
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5. Concluding remarks 

By adopting an appropriate combinatorial perspective, we are able to unify the analysis 
of the number theory and random matrix theory expansions. We reduce showing agreement 
of the two expressions of the n-level density to a combinatorial identity, which we can 
verify for n < 7. As there should be a purely combinatorial proof of this identity, we 
conclude with a few thoughts related to it; we welcome any correspondence with people 
interested in extending these arguments. 

5.1. Verifying the identity formally. We can view the Fourier Identity as a formal iden- 
tity: the indicator functions xa generate a subring C{n) C £oo(]R") that is a quotient of a 
polynomial ring in 2" variables, 

C[xA:AC{l,...,n}]-^C{n), xa^Xa (5.1) 

As a ring of functions, C{n) is certainly reduced, so it is sufficient to check that the identity 
holds over every quotient C{n)/P, where P E Spec C{n) is a prime ideal. 
By equation (14.291) of Lemma l4.6[ this map factors through the quotients 

C'{n) = ^[^A:AC{l,...,n}] ^ ^ Cjn) ^^^^^ 

{xaxb - Xa ■■ A C B) ' J2i^Ai- ■ -xaJ' 

where the second quotient is by the monomials xai ■ ■ ■ xa^. such that xAi ■ ■ ■ Xa^ is identi- 
cally zero as an indicator function in the supported region, namely 

{mi > 0, ...,u„ > 0,^Mi < 2} C M". (5.3) 

(We remark that condition (14.301) does not describe all such products.) 

For C'{n), prime ideals are in one-to-one correspondence with antichains: if {VTi, . . . , Wk] 
is an antichain, the corresponding prime ideal is 

P = {xa- for each i, A^ Wi) + (xa - 1 : for some i,A^ Wi). (5.4) 

Passing to C"{n) just removes 'identically-zero' antichains from consideration. For each 
of the remaining antichains W = {Wi, . . . , W^}, we consider the Fourier Identity under 
the map C{n) — )■ C that sends 

{1 AD Wi for some i, 
(5.5) 
otherwise. 

Verifying that the Fourier identity holds under each of these maps is sufficient to verify the 
full Fourier identity. Assuming C"{n) = C{n) (that is, assuming there are no additional 
relations between the xa)^ this is also a necessary condition. 

We express both sides of the Fourier identity in terms of Euler characteristics. Fix an 
antichain W = {IVi, . . . , 1^^} and let 

k 

S = S{W) = {AC {l,...,n} : ACW;;foTsomei} = [j[0,W;;]. (5.6) 

1=1 
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Thus S is a simplicial set; its vertices are 1J*L^ Wi and its maximal faces are the W^. By 
evaluating as in (15.51) . the sets A with E S evaluate to 1, so the right-hand side of the 
Fourier identity becomes 

E (-l)'^""' = E(-l)""'^'^' = i-ir-'xEAS), (5.7) 

A-.A'^GS AgS 

up to a sign, the Euler characteristic of the simplicial complex S. 

We now express the left-hand side in a related way. First, we determine the value of 
Xaaj + Xa'^aj under the evaluation map. Given J,W C {1, . . . , n}, let U (J, W) be the 
union of segments from V{1, . . . ,n), 

u{j, w) = [w - J, Wuiw- J)] u[wnj,w''u{wn J)] . (5.8) 

It is easy to see the following: 

(1) the two segments are disjoint ifW^0, 

(2) AAJ D W if and only if A is in the first segment, 

(3) A'' A J D W if and only if A is in the second segment. 

In particular, we conclude that, evaluated at W, 

Xaaj + Xa'^aj — S ^ ^, . O-^) 
I otherwise. 

We are only interested in chains where Ai = {1}, so let 

k 

U{J;W) = [jU{J,W,), (5.10) 

i=l 

U{J-W) = nf/(J;W). (5.11) 

Evaluating the inner summand of the Fourier Identity for J gives a sum over all chains in 
U that begin with Ai = {!}. If {1} ^ U, there are no such chains in U, so the sum is 0; 
otherwise, such chains are in bijection with all chains A' in U (with the length off by 1), so 
that 

J2 (-1)"-I^l = (_l)n J2 (-l)l^'l-^ = i-lTxEnm, (5.12) 

chains A chains A' 

Ai = {1}, A,&U 

Aieu 

the Euler characteristic of the order complex Aord(f^)- 

In other words, the Fourier Identity now reads, with the (—1)" canceled and xeuI denot- 
ing the Euler characteristic, 

E -XEnl{U{J-W)) = XEul(5(W)). (5.13) 

JC{2,...,n} 
{1}&U{J;W) 
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Appendix A. Which products of the indicator functions xa are 0? 

In addition to the approach using Euler characteristics in section [STl the authors are in- 
terested in suggestions or answers to the problem of determining which products xa^ ' ' ' Xa^ 
are identically zero in the integration region, where Aj C {1, . . . , n} and 

Xa ■■= x{^^i~ 

ieA ieA'' 

and x{x) is the indicator function of (1, oo). 

In other words, we wish to solve the following linear program: let M be an A; x n matrix 
with each entry ±1, and let 

6 = (l■■•l)^eM^ (A.i) 

c = (1 ■ ■ ■ 1)^ G M". (A.2) 

Minimize c^x = Xj, subject to 

Mx > b, (A.3) 

X > 0. (A.4) 

For j = 1, . . . ,k, let Aj C {1, . . . , n} be the set of +l's in the j-th row of M. Then the 
product xai ■ ■ ■ XAk is identically zero iff one of the following holds: 

(1) the minimum of c^x is 2 or greater, or 

(2) the problem is infeasible. 

The product is nonzero iff the minimum c* G [0, 2). Note that the problem cannot be 
unbounded since c^x > 0. Of course, we could replace the objective function by the 
inequality ^ ■ Xi < 2. 



Appendix B. Sums over chains 

We include a proof of Lemma 14^41 involving sums over chains. The authors thank B. 
UUery for the proof of 6.1(1). 

Lemma B.l. Given chains A, B, we write B ^ AifB is a subchain of A (we include the 
'empty chain' with no sets). If A = Ai d ■ ■ ■ C. Ak, we write k = \A\. Then 

(I) For any n, ^ (^l)'"^' = 

A-.l&Ai, 
Ak = {1, .■■,"} 

i2)Foranyn, ^ (-1)1^1 = (-1)"^^ 

A:Ai = 0, 
Ak = {l,...,n} 

(3) For fixed B with 1 e Bi, ^ (-1)1-^1 = 

AyB,ieAi, 

Ak = {l,...,n} 





ifBl = {l}, 

Otherwise. 
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(1) For n — 1 there is only one possible chain of the desired form, namely {1}. Other- 
wise, there is a bijection between chains A of the desired form with Ai = {1} and those 
with Ai D {1}, by deleting or prepending {1} from the beginning of the chain. Since this 
reverses the parity of \ A\, the sum vanishes. 

(2) Inductively, consider a chain A' : = A[ C ■ ■ ■ C A'^ = { 1 , . . . , n — 1 } on 
{1, . . . ,n — l}. There are t — 1 ways of inserting the element n into the chain while keeping 
Ai = and the last set equals {1, . . . , n}: we can add it into one of the A^, i = 2, . . . ,t, 
or we can insert it immediately after A^ as A^ U {n}, for i = 1, . . . , t. 

The chains A obtained this way contribute (t-l)- (-l)l^'l +t- = 
giving the recurrence 

E (-!)'■" = - E (-1)'^''- (B-i) 

A:Ai=0, A:Ai=0, 
Ak={l,...,n} Afc = {l,...,n-l} 

For n = 1, there is only one such chain, namely C {1}, which has length 2. 

(3) Write B — Bi c ■ ■ ■ C Bk. Choosing Ah Bis the same as choosing k + 1 chains, 
namely, a chain with 1 G Ai and A^ = Bi \ then, for each 2<i</c — 1, a chain from Bi 
to -Bj+i, and a chain from B^ to {1, ... , n}. Thus, we factor our sum as 

E (-1)1^1 = (-1)1^1 ( E E E (-i)'-")' 

^bH.leAi, A:leAi A:Ai=Bi A:Ai=Bk 

Ak = {l,...,n} Ak=Bi Ak=B2 A^={l,...,n} 

(B.2) 

Each Bi is double-counted in the lengths of the chains, so we multiply by (—1)^ = (—1) 
By parts (1) and (2) above, this gives 




^1 ~ ^^H . f_lV-B2|-|Bl|-l...f_nn-|Sfe|-l 



(B.3) 



which is (—1)" when Bi — {1} and otherwise, as desired. 
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